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Abstract 

Transverse-momentum dependent (TMD) light-cone wave functions of a light meson are im- 
portant ingredients in the TMD QCD factorization of exclusive processes. This factorization 
allows one conveniently resum Sudakov logarithms appearing in collinear factorization. The 
TMD light-cone wave functions are not simply related to the standard light-cone wave func- 
tions in collinear factorization by integrating them over the transverse momentum. We explore 
relations between TMD light-cone wave functions and those in the collinear factorization. Two 
factorized relations can be found. One is helpful for constructing models for TMD light-cone 
wave functions, and the other can be used for resummation. These relations will be useful to 
establish a link between two types of factorization. 
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In the collinear QCD factorization for exclusive processes^ such as the form factor of n, 
the nonperturbative effects are included in various light-cone wave functions. At leading twist, 
transverse momenta of partons entering a hard scattering are neglected. The hard scattering can 
be studied in perturbative QCD. In this approach the perturbative expansion for the hard scattering 
often has large corrections, in the form of large Sudakov double logarithms, around the end-point 
regions in which one parton in a hadron carries almost all the momentum of the hadron. These 
large corrections will make perturbative expansions diverge and a resummation is needed for them. 

The solution for resumming these corrections is suggested in [3JII] by taking transverse momenta 
of partons into account, in which one introduces transverse momentum dependent (TMD) light- 
cone wave functions similar to the light-cone wave functions in the collinear factorization. We 
will call the latter as the standard light-cone wave functions. Before a detailed discussion of these 
wave functions some explanation is needed for the nomenclature. The light-cone wave functions 
in the collinear factorization here are called quark distribution amplitudes in the pioneer work 
in pQ. The transverse momentum dependent light-cone wave functions were also introduced in 
the light-front Hamiltonian formulation of QCD and called as wave functions(see the review 0). 
However the TMD light-cone wave functions defined in this letter are slightly different than those 
wave functions, the difference is because of light singularities and will be discussed later. With the 
TMD light-cone wave functions one can make a TMD factorization instead of the collinear one and 
show that the Sudakov logarithms can be resummed. Because the TMD light-cone wave functions 
are generalizations of the standard ones, one may expect to obtain the standard light-cone wave 
functions from the TMD ones simply by integrating the transverse momenta. However, this is 
not the case. In this letter we explore in detail the relationship between these two types of wave 
functions. 

Similar problems also appear in the collinear factorization for inclusive processes like semi- 
inclusive DIS and Drell-Yan process as well as exclusive -B-decays, where large corrections appear 
around edges of kinematical regions. In inclusive processes like Drell-Yan and semi DIS one can 
introduce TMD parton distributions and a TMD factorization for the differential cross-sections 
can be obtained H3 El IE]- in this formalism one finds that the resummation can be conveniently 
performed and there is a factorized relation between TMD parton distributions and the distributions 
appearing in the collinear factorization. In exclusive B-decays, where 6-quark is described by heavy 
quark effective theory, the TMD light-cone wave function of £>-meson can also be consistently 
defined [SJ an d its relation to the light-cone wave function in the collinear factorization is explored 
in With the TMD light-cone wave function one can show that the TMD factorization for 
radiative leptonic decay of S-meson can be verified at one-loop level and the large correction 
around the end-point region can be resummed 11 . Using wave functions or TMD light-cone wave 
functions the form factor of B —* ir transition can also be formulated in a factorized form|12j. 
In other .B-decays it has been shown that the Sudakov logarithms can be resummed (see e.g., 
jl.Sj and references therein). It should be pointed out that TMD factorization is not only helpful 
in resumming large logarithmic corrections but also in probing the 3-dimensional structure of 
hadrons. Various physical quantities can be represented with these wave functions, e.g., various 
form factors [Tl]. generalized parton distributions |15j . single spin asymmetries j!6j . etc. 

TMD light-cone wave functions can be defined from the matrix elements of quark and gluon 
operators in QCD. They have been classified in terms of partonic configurations for various hadrons 
[T7j . In the letter we will take tt as an example, although it is straightforward to extend our results 
to other hadrons. We will use the light-cone coordinate system, in which a vector is expressed 
as a» = (a + ,ar,a±) = ((a + a 3 )/^, (a - a?) / V%, a 1 , a 2 ) and a\ = (a 1 ) 2 + (a 2 ) 2 . We take 
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7r with the momentum P^ = (P + , P~ ,0,0) with P + as the large component and introduce a 
vector u^ 1 = (u + ,u~,0,0). The TMD light-cone wave function of tt can be defined as in the limit 

u + « u~: 

<Mx,i ± ,Cf<) = /^^e"= + ^-^- S M0|9(0)4(co,0) 7 +75i„(~,2)9(2)k(P))|» += „, 

where is the light-quark field. We do not specify if the light quark q is a u- or <i-quark and 7r 
is 7T° or 7r + , which is not important for our discussion. L u is the gauge link in the direction u: 

L u (oo, z) = Pexp ^— ig s J dXu ■ G(Xu + z)^j . (2) 

It should be understood that the contributions proportional to any positive power of u + /u~ are 
neglected. This definition was first proposed in [3]. It is gauge invariant in any non-singular gauge 
in which the gauge field is zero at infinite space-time. The TMD wave function has an extra variable 
( 2 beside the momentum faction x, the transverse momentum k± and the renormalization scale 
ft. The evolution of this variable will generate the resummed Sudakov logarithms as shown in [3]. 
This will be confirmed in this letter. The evolution with the renormalization scale /i is simple: 

d(t> + (x,k ± ,(,fi) 

V q- = 2-y q (f) + (x,k ± , (3) 

where -y q is the anomalous dimension of the light quark field q in the axial gauge u ■ G = 0. In 
the axial gauge, the gauge links in our definition disappear. Some general features of TMD light- 
cone wave functions defined by using the light-cone gauge link with u + = have been studied 
in |18l I19j . In |18j the operator product expansion was employed to analyze the short distance 
behavior of wave functions. However, the TMD light-cone wave functions defined with u + = will 
have light-cone singularities like 1/(1 — x), as pointed out in j^H] for TMD parton distributions. 
We will also show through our calculation that these singularities exist if one sets u + = at the 
beginning. With a finite, but large the light-cone singularities are regularized. Because of the 
small but finite u + , the TMD light-cone wave function is not the wave function introduced in the 
light-front Hamiltonian formulation of QCD[Sj. However, with our definition one can still derive the 
so-called Drell-Yan-West relation(See ESI)) which shows that the behavior of the TMD light-cone 
wave function of a hadron in the region of x — > 1 is related to the structure function of the hadron 
in the region of x — > 1 and the leading power behavior of the form factor factor of the hadron. 
With the definition of TMD light-cone wave function it has already been shown that the TMD 
factorization for the form factor in the transition of 7*7r — > 7 can be consistently factorized at least 
at one- loop level |211. It should be noted that the TMD light-cone wave function defined in Eq.(l) 
is real. This can be shown by using parity- and time-reversal transformation and with the fact that 
the TMD light-cone wave function depends on the vector u through 
The standard light-cone wave function is defined aspQ: 

<Mx, M ) = / ^ r e jfc+2 ~(0|g(0)L] l (oo,0)7 + 75^n(oo,z)g(z)|vr(P))U +=0iZ -. ±=0 , (4) 

where the gauge link is along the light-cone direction = (0, 1, 0, 0). Comparing the two definitions 
one would expect in the limit u + — > that the standard light-cone-wave function can be obtained 
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by: 

$ + (x,fj,) = J d 2 kx(f) + (x,k±,C,fj,), (5) 

where the integration over k\ is from to do. However, this is not true. The reason is as follows: 
By the generalized power counting rule [22], <fi+ is proportional to fcj as k\ — > oo. Hence the 
transverse momentum integral is ultraviolet (U.V.) divergent. In Eq.(4) the integration over the 
transverse momentum is supplemented with systematic U.V. subtractions and this generates a 
renormalization scale dependence in <3?_|_. This leads to the Efremov-Radyushkin-Brodsky-Lepage 
evolution equation [2D- In the above integration one may also implement an U.V. subtraction by 
introducing a suitable cut-off for k±. This is easy to do at one- loop level, and but difficult to extend 
beyond one-loop. Also the limit u + — > or £ — > oo is nontrivial. 

Although it is tricky to establish the above relation, the two types of wave functions are related 
to each other in other ways. If we transform the TMD light-cone wave function into the impact 
space 

<t> + (x,b,(,») = J d 2 fc±e^0+(a;,fcj.,C,/i) (6) 
it can be shown that a factorized relation exists between two wave functions if b is small: 

0+(x,6,C,/i) = / dyC b (x,y,C,b,ii)$+(y,n)+0(b), (7) 
J o 

where the function C can be calculated in perturbative QCD and does not have any soft divergence. 
When b is small, it corresponds to a large momentum scale, and hence its dependence must be 
calculable in perturbation theory. Q can be understood as a large scale, its dependence shall also 
be perturbative. The factorization theorem asserts that all nonperturbative effect in (/)+(x, b, CjM) 
is resided in &+(x, fx) for small b. We will show this is true at one-loop level and our result can be 
extended beyond one-loop. 

Another interesting relation can be found if k± is large, which can be generated by exchanges 
of gluons between partons inside tt and these gluons are hard. Hence the behavior at large k± can 
be studied with perturbative QCD. One expects the type of factorization: 

1 f 1 

<t>+(x,k L ,(,v) = j o dyC±(x,y,C,n)$+(y,n) + 0(/cJ 4 ). (8) 

The factor kj 2 is determined by the power counting rule in |22| . 

As discussed above, the two functions C b and C± can be calculated in perturbative QCD. To 
do that, we take a partonic state to calculate the two wave functions at one-loop order, in which we 
have both infrared- and collinear divergences. We regularize infrared divergences by taking a small 
gluon mass A, and collinear divergences by taking a small quark mass m q . If the factorization in 
Eq.(7) is correct, C b will have no singular dependence on these small masses. It should be noted 
that our results for the two wave functions at one-loop level will also be useful to establish TMD- 
and collinear factorization theorems for exclusive processes at one-loop. We take the partonic state 
\q(k g ), q{kq)) to replace ir in the above definitions, the parton momenta are given as 

k£ = (k+,k-,k q± ), A# = (k+,k=,-k q± ), k+ = x P + , k+ = (l-x )P + = x P + . (9) 

These partons are on-shell. At tree-level the wave functions are trivial: 

(f> { +\x,k ± X) = 5(x-x )5 2 (k ± - k q± )fa, < 0) (x )/U ) = 5(x-x o )0o, (10) 
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where <f>Q is a product of spinors 4>o = v(kq)^ + ^^u(k q ) / P + . We will always write a quantity A 
as A = A^ + A^ + • • •, where A^ and A^ stand for tree- level- and one- loop contribution 
respectively. 

At one- loop level, one can divide the corrections into a real part and a virtual part. The real 
part comes from contributions of Feynman diagrams given in Fig.l. The virtual part comes from 
contributions of Feynman diagrams given in Fig. 2, which are proportional to the tree-level result. 




Figure 1: The real part of one- loop contribution to the TMD light-cone wave function. The double 
lines represent the gauge links. 

The contribution from Fig. lb and Fig.lc and Fig. Id to (/>+ reads: 
, , , ... 2a s f„, . / 1 \ x 1 , 1 . q 2 \ + \ 2 \ , 

X 

A q = (1 - y) 2 ™ 2 - + 2A 2 + ( k ± - yk q ±) 2 , y = —, q± = k ± - k q± , 

Xq 

4>+(x,k ± ,()\ic = {^(x - x ) (z—Zr) n 1 I \ 6 ( x ~ x o)zT~Tv2 ln 9± + A A ~ 

2 X 



.x-x J+(l-x Q )A q 2^ VJ qi + \i ( 2 4 
&q = (1 - y) 2 m\ + y\ 2 + (k ± - yk q± ) 2 , y = 



x 
2a 1 

<j)+(x,k ± ,C)\ld = ~^f A 2 + g 2 5 ( X ~ X o)4>0, (11) 

where we have already taken the limit £ — > oo and only kept the leading terms. If we set u + = or 
C = oo at the beginning, the contribution from Fig. lb and Fig.lc will be divergent as (x— xo) -1 when 
x — > xq, and Fig. Id will not lead nonzero contribution. The divergence is the light-cone singularity 
mentioned before, caused by the exchanged gluon if it has the momentum = (kg ,k~ ,k g ±) 
with a vanishing small kg. With a finite, but large £ the divergence is regularized. It results 
in the +-distributions and the terms with ln£ in the above expressions. From the above results 
one clearly sees that there are collinear- and infrared singularities when transverse momenta are 
small. It should be noted that the contributions from Fig. lb and Fig.lc are related to each other 
by charge-conjugation. 

The contribution from Fig. la is complicated. But for the function C± and Cb we only need the 
leading part of the contribution with k± — > oo. The leading part can be written as: which is 

2d I i x 1 — x \ 

M^k ± ,Q\ la = ^ 1 T^\-0( Xo -x) + —0{x-x ))<h 

+" finite terms" , (12) 

where we introduce a "cut-off" A, which is also contained in the second line, so that the total does 
not depend on it. The second line behaves like fej 4 with k± — > oo. One can show that the second 
line is irrelevant here, but it is relevant to what is neglected in Eq.(7,8). 
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Figure 2: The virtual part of one-loop contribution to the TMD light-cone wave function. The 
double lines represent the gauge links. 



The virtual part of the one-loop correction is from the Feynman diagrams given in Fig. 2. 
Contributions from each diagrams are: 



(f>+(x,kj_,C)\ 2a 

4>+(x,k ± ,()\2b 
4>+(x,k±,()\2c 



cj) + (x,k±,()\2d = </)+ > (x,k±,C) ■ ^ 



ln-^ + 21n-^ 



mi 



(f) + (x,k±,()\2e = (f>+\x,k±,() ~ln^, 



07T 



2 In At + 2 In 



c 2 4 



mi 



mi 



m„ 



2^2 



+ In 2 -± - In 

A z 



A 2 



2tt 2 



+ 4 



4> + (x,k ± ,()\2c\x ^x - 



(13) 



These results also contain collinear- and infrared divergences. The complete one-loop contribution 
<fr+ is the sum of contributions from the 10 Feynman diagrams in Fig.l. and Fig. 2. From the above 
one-loop result one can already obtain the relation in Eq.(8) with the function C± determined as: 



2a s 
3vr 2 



-9{y-x) 



x 1 



yx-yj + 



+ 0(x-y) 



+6(x - y) 



k 



In -± + ln(l - x) 2 + In x 2 - 1 



1-x 1 
1-yx-y 

+ 0(a 2 s ). 



+ 



1-2/ 



(14) 



In the above the +- prescription acts on the distribution variable y. We have used the identity 



x 9{y — x) 
y x-y 



x 9(y — x) 



V x-y J + 



5{x-y) [l + ln(l-x)], 



(15) 



where the +- prescription in the left hand side acts on the distribution variable x, while the +- 
prescription in the right hand side acts on the distribution variable y, 

To derive the relation in Eq.(7) one needs to transform the above results into the impact 
parameter space and to obtain one-loop results of <I> + . The Fourier transformation can be done 
straightforwardly. We find that there is a cancelation of infrared divergence between the real- 
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and virtual part correspondingly. We present our results in the combination in which infrared 
divergences are canceled: 



cp + (x,b,C,fj)\ 



2a A 



16+2c 



1 



——4>o < -5(x - x ) 



3vr 



+9(x - x) 



In 

2™2 



Cx z ) + 21n(6V) + 21n(^C^o) - 4 - 7r 



2/-2„.2\ 



x ln(6 2 m 2 (l - y) 2 ) 



:r 



+ 0(6), 



4> + (x,b,C,fJ-)\ 



2a.s 



lc+2/ 



3tt 



1 



— </> ^ -<5(x - x ) 



X — Xq 

- ln 2 (b 2 ( 2 x 2 ) + 2 ln(&V) + 2 ln(& 2 C 2 * 2 ) - 4 - vr 5 



7 (a; - x ) 



1-x ln(6 2 m 2 (l-y) 2 ) 



1 - x 



X — Xq 



+ 0(b), 



(j)+(x,b,(,fJ-)\ld+2b+ 



2c 



^<j>o6(x-x Q )]n(b 2 fi 2 ) + 0(b), 



(16) 



with b = be 7 /2. The contribution from Fig. 2a and Fig. 2d in the b-space can be easily read off from 
Eq.(15) With these results in the 6-space one can derive the evolution of £: 



d 

C^4>+{x,b,C,fi) 



2a Sl CWe 2 ?" 1 2a s , Q 2 (l - xftfe 2 ^ 1 

in m 

3vr 4 3vr 4 



fa* 
" 3vr 



ln(&V) + In 



C 2 x(l — x) 
e/i 2 



<j) + (x,b,(,v) + 0(a 2 ) 
<P + (x,b,(,v) + 0(a 2 ). (17) 



This equation was derived first in Our result agrees with that in 3 . The solution of the 
equation will contain the Sudakov logarithms in a resummed form. 

The one- loop contributions to <£+ are represented by the same diagrams given in Fig.l and 
Fig. 2, except those in which the gluon is exchanged between gauge links. Some of them contains 
light-cone singularities because the gauge link is along the direction re, i.e., u + = is set at the 
beginning. This singularity is canceled between the real- and virtual part correspondingly. We 
present our results in the combination free from the singularity: 



*+(x,/i)|i c + * + (s,/i)| 2 / 



2a f 



3tt 



<t> Q 9(x - x ) 



1 — x 



2a, 



$+0,^)116 + $+(x,[i)\ 2c = —4>o9(x -x) 

07T 



(1 - X )(X - Xq) 
X 



m 2 q (l - y) 2 



X (x - X ) 



m 2 {l-yf 

m — — 5 

u 2 



(18) 



Again, the contribution to from Fig. la is complicated. However the relevant part is just by 
integrating the leading part of <^>+|i a over k± with dimensional regularization. We have the relevant 
part as: 



2a s 
3vr 



hi 



A 2 



1 



— 6(x - X) + — — 6{X - Xq) 

Xq I-Xq 



+" finite terms" , 

where the factor —1 comes from manipulation of 7-matrices in d-dimension. 



(19) 
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With the above results we can extract the function Cb(x,y,(,b, (J,). At tree-level 



5(x - y). 



(20) 



At one-loop we have: 



2a, 
3vr 

+-6(x-y) 



(ln(&V) 



x 1 
-6(y-x) + - 

y i 



-6{x - y) 



In 



2^2 a2„.2\ 



CV)-lnW(l-y) 2 )-8-2vr 



2(ln y 2 + ln(l - yf) ln(6 2 /i 2 ) + 2 hi(6 2 C V) + 2 ln(6\ 2 (l - yf) 



2/-22\ 



(12 a2/ 



+0(2/ -X) 



x ln(6 2 /i 2 ^ 



y sc-y 



(x - y) 



l-x\n(b 2 ^ 
1 - y x-y 



(21) 



In the above the +-prescription acts on the variable y. Clearly, it is free from any collinear- or 
infrared singularity as expected. 

To summarize: In this letter we have performed a study of the TMD light-cone wave function 
of a 7r meson which appears in TMD factorization of exclusive processes. We have established two 
factorized relations between the TMD- and the standard light-cone wave function, the latter is 
relevant in collinear factorization of exclusive processes. One relation is that the TMD light-cone 
wave function can be written as a convolution of the standard one with a perturbative coefficient 
function Cj_ when the transverse momentum is large. This relation is helpful for constructing models 
of the TMD light-cone wave function. Another is that the TMD light-cone wave function in the 
impact parameter b space can be written as a convolution of the standard one with a perturbative 
coefficient function Cb when b is small. The function Cb contains only perturbative effect and is 
determined at one-loop level. This factorized relation can be extended beyond one-loop level and 
it is useful for resummation of Sudakov logarithms. From our results we confirm the result of the 
evolution equation of £, first derived in [3]. The solution of the equation resums the large Sudakov 
logarithms. 
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